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ABSTRACT 

We have examined a set of pure luminosity evolution (PLE) models in order to explore 
up to what extent the rapidly increasing observational constraints from faint galaxy 
samples can be understood in this simple framework. We find that a PLE model, 
in which galaxies evolve mildly in time even in the rest frame UV, can reproduce 
most of the observed properties of faint galaxies assuming an open (f2 ~ 0) universe. 
In particular, such a model is able to fit reasonably well the number counts in the 
U, bj, rf, I, and K bands, as well as the colour and redshift distributions derived 
from most of the existing samples. The most significant discrepancy between the 
predictions of this model and the data is the z distribution of faint if-selected galaxies. 
Significantly worse fits are obtained with PLE models for the theoretically attractive 
value of r2 1 , although a simple number luminosity evolution model with a significant 
amount of merger events fits the data also in this cosmology. 

Key words: galaxies: evolution - galaxies: photometry - galaxies: redshifts - cos- 
mology: miscellaneous. 



1 INTRODUCTION 

An impressively large amount of observations of faint field 
galaxies has been collected over the last fifteen years. The 
determination of galaxy number counts, colour and redshift 
(z) distributions of increasingly deep samples is now part 
of the observational routine of several established groups. 
Since the review paper by Koo & Kron (1992), galaxy pho- 
tometric surveys have been extended to bj = 27.5 (Metcalfe 
et al. 1995a) and to A' < 24 (Gardner, Cowie, & Wainscoat 
1993; Soifer et al. 1994; Djorgovski et al. 1995), and spec- 
troscopic surveys to A" < 20 (Songaila et al. 1994), / < 22 
(Crampton et al. 1995) and to B < 24 (Glazebrook et al. 
1995a). These data provide useful constraints for both cos- 
mological models and galaxy evolutionary models. 

As soon as deep optical counts of galaxies became avail- 
able, it was realized that no-evolution (uE) models, i.e. 
models in which the absolute brightness and the spectra 
of galaxies do not change in time, predict a surface density 
of galaxies at faint magnitude significantly lower than the 
observations (see the review by Koo & Kron 1992). The ex- 
cess in the observed number counts with respect to the nE 
model predictions is of the order of ^ 4 to 5 at bj ~ 24 and 
~ 5 to 10 at bj ~ 26 (Maddox et al. 1990; Guiderdoni & 
Rocca-Volmerange 1991, hereafter GRV91). 

Simple pure luminosity evolution (PLE) models, which 



allow for brightness and spectral evolution in the galaxy pop- 
ulation, were shown to provide a better fit to the faint galaxy 
number counts (Tinsley 1980; Bruzual & Kron 1980 (here- 
after BK80); Koo 1981,1985). However, these early PLE 
models were ruled out on the basis of comparisons with the 
results of z surveys of faint galaxies (Broadhurst, Ellis, & 
Shanks 1988; CoUess et al. 1990; Koo & Kron 1992 and 
references therein), which failed to reveal the large number 
of high z galaxies predicted by the models. The z distribu- 
tions in these surveys appeared to be in agreement with the 
predictions of nE models. 

In an attempt to explain all aspects of the data with 
a single model, this apparent paradox (counts require PLE 
model but z distribution only admits nE model) prompted 
the development of a number of less straightforward mod- 
els. These models relax some of the assumptions adopted in 
nE and PLE models, such as, the constancy of the volume 
number density of galaxies, and/or the standard cosmology. 
The excess observed in the faint counts can be explained, 
for example, by either non-conservation of the number of 
galaxies due to merger events (GRV91; Broadhurst, Ellis, & 
Glazebrook 1992; Carlberg & Chariot 1992) or dwarf galax- 
ies which have faded and/or disappeared in recent epochs 
(Broadhurst et al. 1988; Cowie, Songaila & Hu 1991; Babul 
& Rees 1992). Fukugita et al. (1990) have proposed to revise 
the cosmological model introducing a non-zero cosmological 
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constant. 

More recent and refined versions of PLE models have 
shown that at least some of the discrepancies with existing 
data can be substantially reduced also in the framework of 
these models (Guiderdoni & Rocca-Volmerange 1990 (here- 
after GRV90); Gronwall & Koo 1995 (hereafter GK95); Met- 
calfe et al. 1991,1995a). 

The main goal of this paper is to explore up to what 
extent the properties of the faint galaxy samples can be un- 
derstood in the framework of PLE models, without invok- 
ing the more exotic scenarios proposed in the literature. We 
examine a new set of PLE models and compare their pre- 
dictions with the most recent data, including counts in all 
available photometric bands, and colours and z distributions 
as a function of magnitude for various samples, selected in 
different bands. A global comparison of this kind is essen- 
tial because each different data set provides constraints on 
different PLE model parameters. 

We find that standard PLE models built along the lines 
described by BK80, but with mild gala^cy evolution, similar 
to that assumedhy Metcalfe et al. (1991, 1995a), but derived 
here consistently from the Bruzual & Chariot (1993, here- 
after BC93) spectral evolution models, provide satisfactory 
fits to most of the observational data in an O ~ universe. 
The only data set which our models do not reproduce well 
is the z distribution of the it'-selected sample for K > 17 
by Songaila et al. (1994). Despite this failure, and before 
new samples either confirm or modify the Songaila et al. re- 
sults, we think that it is useful to compare the faint galaxy 
surveys with the predictions of simple PLE models in or- 
der to constrain the range of evolutionary and cosmological 
parameters used in these models. 

In §2 we describe the main ingredients of our PLE 
model, i.e. the shape and normalization of the luminosity 
function of galaxies of various morphological types and the 
adopted models for their spectral evolution. In §3 we com- 
pare the predictions of our models with the observational 
data. The conclusions of our work are presented in §4. 

2 STANDARD PLE MODELS 

Standard PLE models (BK80; Metcalfe et al. 1991) assume 
that galaxies maintain at every z the same proportion of 
the various morphological types as observed locally (2 = 0), 
and that the spectral evolution of these galaxies is well de- 
scribed by models which at low z reproduce the colours and 
^-corrections determined from observed galaxy spectra. As- 
suming a geometry for the universe, and scaling the galaxy 
luminosity function (LF) to match the number counts at a 
given magnitude in a specific band, one can compute the 
expected galaxy number counts in various bands, as well as 
the colour and z distributions as a function of magnitude. 

2.1 Modelling faint galELxy counts and colour and 
z distributions 

In a homogeneous and isotropic universe, the number of 

galaxies of each type brighter than a given magnitude m 
can be calculated from the integral 

Ni{<m)= / / ct>i{M)dM^^-^j^dz, (1) 

Jo J M,nin 



where 4>i{M) and dV{z,Q) are, respectively, the local LF 
for the type and the comoving volume element. In this pa- 
per we will consider only the standard (A = 0) Priedmann 
cosmology defined by Ho and Q. We use Hq = 50 km s~^ 
Mpc^^ throughout this paper. 

The total number of galaxies brighter than m, N{< m), 
is obtained by adding A^i(< m) over all types considered. 
From (1) we obtain the z-distribution Ni{< m;z,z + dz) 
by integrating over the specific z range {z,z + dz). In (1) 
Zmax = min (zf,ZL), where z/ is the assumed z of galaxy 
formation and zl is the value of z at which the Lyman 
continuum break is shifted to the effective wavelength of 
the filter being considered, and the galaxy presumably be- 
comes dark (Madau 1995). We use the following filters: U 
(Koo 1981), B and V (Buser 1978), bj and r/ (Couch & 
Newell 1980), / and K (Wainscoat & Cowie 1992). For the 
U, bj, rf, I, and K bands, zl — 3, 4, 6, 8, and 23, respec- 
tively. 

The integration over dM in (1) extends up to 
Mmax{m, z)=m- 5/og '^^^^'j^°'^^ - corr{z\ Ho, fi, Zf),{2) 

where dL{z;Ho,0.) is the luminosity distance measured in 
pc, and corr{z; Ho, n, Zf) is the correction needed to obtain 
the galaxy rest frame magnitude from its observer frame 
magnitude. In the riE model, it is just the fc-correction. In 
the PLE model, it is given by the (e -|- fc)-correction, which 
also includes the effects due to the intrinsic galaxy luminos- 
ity evolution. We computed this correction up to 2 = Zmax 
from the synthetic spectral energy distribution (SED) of the 
various galaxy types for each photometric band. 

Finally, the colour distribution Ni{< m; c, c+dc) can be 
derived from the z distribution Ni{< m;z,z + dz) for each 
type of galaxy by using the cosmology dependent relation 
between colour and z, Ci{z; Hq, Q, Zf), given by the adopted 
spectral evolution model. The colour distributions for the 
different types are then added up. 

2.2 Luminosity function 

There is increasing evidence that the galaxy LF varies with 
galaxy type. Bingelli, Sandage & Tamman (1988) have de- 
rived morphology dependent LFs for galaxies in the Virgo 
cluster. Efstathiou, ElUs & Peterson (1988), Shanks (1990) 
and Loveday et al. (1992) have shown that the LF in 
the field is well described by the analytical expression of 
Schechter (1976), with values of a and M* which depend 
on the galaxy morphological type. In general, blue galaxies 
show a steeper slope than red galaxies. 

We adopt the values of a and Mj* derived separately 
for early and late type galaxies by Efstathiou et al. (1988) 
from the Anglo- Australian Redshift Survey data, listed in 
Table 1. We assume that these parameters are valid from 
Ml,. = -24 to Ml,. = -15.5 (cf. Loveday et al. 1992). To 
obtain the local LF in the U, rf, I, and K bands we shift 
M^. according to the colours at 2; = for each galaxy type 
listed in Table 1. The K band LF thus obtained is consistent 
with the observed one (Mobasher et al. 1995; Glazebrook et 
al. 1995b). 

We have not used the more recent determination of a 
and . by Loveday et al. (1992), because of the bias men- 
tioned by these authors against identifying early-type galax- 
ies in their sample. As shown by Zucca, Pozzetti & Zamorani 
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Table 1. Local galaxy mix, LF parameters", and local colors 



Type 


Fraction'^ 


a 






B-V 




bj - rf 




bj-K 


E/SO 


0.28 


-0.48 


-20.87 


0.95 


0.95 


0.75 


1.55 


2.32 


4.13 


Sab-Sbc 


0.47 


-1.24 


-21.14 


1.15 


0.68 


0.32 


1.21 


1.83 


3.47 


Scd-Sdm 


0.22 


-1.24 


-21.14 


0.54 


0.43 


-0.05 


0.88 


1.39 


2.96 


very Blue (vB) 


0.03 


-1.24 


-21.14 


0.12 


0.07 


-0.49 


0.32 


0.63 


2.06 



" a, M* from Efstathiou et al. (1988); *| in units of 10"^ Mpc-^ (§2.3). Hq = 50 km s-^ 

Mpc^-*^ throughout this paper; 

From BC93 {tg = 16 Gyr, Zf = 4.5, H = 0); 
= Mix from Ellis (1983). 



(1994), the correction for this bias moves the Loveday et al. 
values of (a, M*) towards those of Efstathiou et al. (1988). 

2.3 Count normalization 

The uncertainties in the normalization of the local galaxy LF 
translate into a major source of uncertainty in the models 
for faint galaxy counts. Moreover, different selection effects 
in photographic and CCD data may lead to difficulties in 
comparing bright and deep surveys (McGaugh 1994). The 
reduced size of well calibrated samples at bright magnitude 
and the local fluctuations caused by galaxy clustering lead 
to large variations in the bright counts. In principle, this 
problem would be avoided by using the APM galaxy counts 
in the range 15 < bj < 20.5 derived from 4300 dog^ of the 
sky (Maddox et al. 1990). However, while these counts 
are close to the mean of previously published data at bj ~ 
19 — 20, they show a steeper slope for bj < 19 than previous 
determinations. Thus, the APM counts at bj < 17 are below 
models normalized to match the counts at 6j = 19 — 20. 
Metcalfe, Fong & Shanks (1995b) show evidence of a possible 
magnitude scale error in the APM galaxy survey, whose size 
could be large enough to cause the apparent disagreement 
between the APM galaxy counts and the predictions from 
standard models in the range 17.5 < 6j < 20. 

Because of this uncertainty in the counts at bright mag- 
nitude, we have chosen to scale our model predictions to the 
observed number counts in the range 19 < bj < 19.5, namely 
log No = 1.98 gal deg"^ (0.5 mag)"^ (sec also GRV90). This 
normalization is consistent with that obtained recently for 
the K band LF (Glazebrook et al. 1995b), as well as for the 
bj band (Ellis et al. 1995). 

Finally, in order to compute $*, i.e. the normalization 
of the LF for each morphological type, we have adopted the 
local galaxy mix derived by Ellis (1983) for bj < 16.5 from 
the DARS data (see column 2 in Table 1). This galaxy 
mix is in good agreement with the percentage of ellipti- 
cal and spiral galaxies derived recently by Efstathiou et al. 
(1988). The same values of $* are used to model different 
photometric bands. The bj LF summed over all types de- 
rived from our models is well described by the parameters 
a = -1.2, Mj. = -21.1, = 2.5 x 10"^ Mpc"-\ 

2.4 Spectral energy distributions 

Our galaxy spectral energy distributions (SEDs) are based 
on the BC93 galaxy spectral evolution library. The BC93 



models axe built from a library of stellar tracks which in- 
cludes all evolutionary stages for stars of solar metallicity. 
Empirical noar-UV to near-IR spectra of galactic stars, ex- 
tended to the far-UV by means of model atmospheres, are 
used in the synthesis. The inclusion in these models of the 
Post-AGB stellar evolutionary phase, which contributes sig- 
nificantly in the UV spectral range in old stellar populations 
(Magris & Bruzual 1993), represents an improvement over 
previous spectral evolution codes. 

We used the following procedure to select a subset of 
galaxy models from the BC93 library. First, we selected 
models on the basis of their ability to reproduce the shape 
of the continuum and spectral features in the SEDs of local 
galaxies of different morphological types, and verified that 
the fe-corrections computed from the models were in agree- 
ment with the empirical ones. Then, among these models 
we selected by trial and error those which reproduced better 
the empirical constraints on galaxy evolutionary properties, 
namely, galaxy number counts and colour and z distribu- 
tions. 

Thus, we have to compute the number count model as 
outlined in §2.1 in order to find which galgixy spectral evolu- 
tion models produce the best fit to the counts and colour and 
2 distributions. The feedback in this procedure is unavoid- 
able due to the lack of independent constraints on galaxy 
evolution (BK80, GK95). 

2.4-1 Evolutionary constraints on the SEDs 

The following aspects of the data guided our choice of the 
SEDs. Strong evolution for 2; < 1 is ruled out for all galaxy 
types by the observed z distributions of Broadhurst et al. 
(1988) and CoUess et al. (1990, 1993), which are close to 
the nE prediction up to bj = 22.5. CoUess et al. (1993) esti- 
mate an upper limit to the amount of luminosity evolution 
of AMb ~ —1.2 for z < 1, providing a strong constraint to 
the evolutionary corrections. The observed number counts, 
particularly in the bj band, require optically mildly-evolving 
SEDs, over the entire z range, in agreement with Koo, Gron- 
waU, & Bruzual (1993, hereafter KGB93). Additional evi- 
dence in favor of mild evolution derives from recent work on 
the z evolution of the LF in B-, I-, and A'-selected samples 
by CoUess (1995), Lilly et al. (1995b), and Glazebrook et 
al. (1995b), respectively. 

We introduce the spectral class of very Blue (vB) galax- 
ies in order to explain the bluest colours, bj — rf ^ 0.3, ob- 
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Table 2. Initial mass function" 



Table 3. Galaxy SED model parameters 



IMF 



mi m2 



Scalo (1986) 



-2.60 
0.01 
1.75 
1.08 
2.50 
1.63 



0.10 
0.18 
0.42 
0.62 
1.18 
3.50 



0.18 
0.42 
0.62 
1.18 
3.50 
125 



100.530 
1.14430 
0.25293 
0.34842 
0.44073 
0.14819 



Salpeter (1955) 1.35 0.10 125 0.17038 

"/<("^) = Cim~(^+'^«) for mi < m < m2, 
mi and m2 in Mq units. 

served in the deepest 6j-selected surveys (Glazebrook et al. 
1995a). vB galaxies are meant to reproduce a population of 
starburst galaxies present at each z, whose evolution does 
not follow the pure luminosity prescription. Wc assume that 
star formation in these galaxies keeps their SED constant in 
time. For vB galaxies we adopt the LF of spiral galaxies. 
The observed fraction of galaxies bluer than B — V = 0.6 at 
the faint magnitude limit 10%), is reproduced assuming 
that vB galaxies represent ^ 3% of the total local mix (Ta- 
ble 1). Thus, at the few % level, our PLE models include a 
fraction of non-passively evolving galaxies. 



3 RESULTS 

3.1 Reference models 

A BC93 galaxy evolution model is specified by the star for- 
mation rate (SFR) ?/>(t) and the initial mass function (IMF). 
Once Ho, Q,, and Zf arc specified, the age tg of the model 
SEDs and their observer frame properties are fixed. We 
considered models with tg = 16 and 12.7 Gyr, correspond- 
ing to Zf = 4.5 and 10 for f2 = and 1, respectively. Wo 
tested models computed for the Salpeter (1955) and the 
Scalo (1986) IMFs (Table 2). Table 3 defines the BC93 
model SEDs selected as described in §2.4 that will be used 
to build our reference models for the counts as indicated in 
§2.1. A brief comment on these SEDs follows. 

Models for E/SO and Sab-Sbc galaxies are characterized 
by the SFR ■tp{t) oc exp(— t/r), where t is the e-folding 
time characterizing this form of ip(t). The spectral prop- 
erties of nearby E/SO galaxies are reproduced well by both 
the T = 1 Gyr model (hereafter ri model) and a model in 
which star formation takes place at a constant rate during 
the first Gyr in the life of the galaxy (hereafter Bi model) 
for either the Salpeter or the Scalo IMF, and by the r = 2 
Gyr model (or T2 model) for the Scalo IMF. The observed 
bj — r-f and B — K colour distributions are reproduced more 
closely by our f2 = reference PLE model if we represent 
the E/SO galaxy SEDs by the Ti and T2 models, rather than 
by the Bi model (Table 3). 

The local properties of Sab-Sbc galaxies are well de- 
scribed by the r = 4 Gyr (or T4) Salpeter IMF model, or 
the r = 10 Gyr (or no) Scalo IMF model. The Scd-Sdm 
galaxies are described satisfactorily by a model in which 
stars form at a constant rate following the Salpeter IMF 
(hereafter cons model). The SED of the vB galaxies can be 
approximated by a model with constant SFR seen at an age 



n 


Type 


SFR 


IMF 


tg (Gyr) 





E/SO 


Ti,r2 


Scalo 


16 




Sab-Sbc 


no 


Scalo 


16 




Scd-Sdm 


cons 


Salpeter 


16 




vB 


cons 


Salpeter 


0.1 


1 


E/SO 


Si,ri 


Scalo 


12.7 




Sab-Sbc 


T8 


Scalo 


12.7 




Scd-Sdm 


cons 


Salpeter 


12.7 




vB 


cons 


Salpeter 


0.1 



E/30 ^ , 




Figure 1. fc-correction in the Tf band for galaxies of different 
morphological types from BC93 models (dotted lines, sc: Scalo 
IMF; sp: Salpeter IMF) and from the observed SEDs described 
in the text (solid lines). 



of 0.1 Gyr after the last major event of star formation. This 
SED is even bluer than that of NGC 4449 (BC93). 



3.1.1 Scalo IMF vs. Salpeter IMF 

The constraints mentioned in §2.4 led us to adopt Scalo IMF 
models for early- type galaxies (Table 3). The Scalo IMF is 
less rich in massive stars than the Salpeter IMF because of 

the steeper slope of the former at the high mass end. The 
high number of massive stars in the Salpeter IMF models 
produces a large amount of UV flux at early times, making 
high z E galaxies detectable in current deep surveys and 
producing an excess in the faint bj counts with respect to 
the observations. An altcruative way to reduce the UV flux 
is to assume a significant amount of dust extinction inside 
the galaxies (GK95). 

There is some ambiguity in the literature on the func- 
tional form of the Scalo IMF. Following a suggestion by D.C. 
Koo, BC93 have divided the IMF given by Scalo in graphical 
form into the 6 different segments listed in Table 2. For the 



-segment the IMF is given by fi(rn) = am 



The 



constants c; are computed from the requirement of continu- 
ity of the IMF (cf. Guiderdoni & Rocca-Volmerange 1987, 
hereafter GRV87). 



Table 4. Intrinsic luminosity evolution 





2300A 


1150A 


11000 A 


5500A 


Type 










E/SO 


-1.6 


-6.0 


-0.6 


-2.3 


Sab-Sbc 


-1.0 


-1.6 


~ 


~ 


Scd-Sdm 


~ 


~ 


-1-0.4 


-1-0.6 



" Rest frame wavelength 



3.2 k and (c + fc)-corrections 

Fig 1 shows the model and the empirical fe-corrections in 
the rf band. These models, listed under 51 = in Table 
3, reproduce the flattening at high z shown by the empiri- 
cal ^-corrections of ellipticals and spirals (cf. Cowie et al. 
1994). In contrast, the fc-corrcctions of ellipticals modeled 
by Rocca-Volmerange & Guiderdoni (1988) continue to in- 
crease with z. This difference is mainly due to the lower UV 
flux in the GRV87 model for E galaxies, resulting from the 
lack of Post-AGB stars in their population synthesis. 

The fc-corrections derived from the synthetic spectra arc 
in agreement with the fc-corrections in the {U, bj,rf,I) bands 
computed from the observed spectra of Pence (1976) up to 
z = (0.6,0.9,2.2,2.5), with the compilation in bj by King 
& Ellis (1984) up to 2 = 1.5, and with the K and B band 
estimates by Cowie ct al. (1994) up to 2 < 3. Our E/SO 
models reproduce also the fe-correction at higher z computed 
from the average observed E galaxy SED of BC93. 

Figures 2a and 2b show the k and (e -I- fc) corrections 
in the bj and K bands, respectively, for the SEDs selected 
in the SI = case (Table 3). We see here the flattening of 
the (e -|- fc)-corrections in bj and r/ at high z, which was 
assumed entirely ad hoc by Metcalfe et al. (1991, 1995a). 
If we use Salpeter IMF models for the Ellipticals, the (e -|- 
fc)-correction in bj continues to decrease with z, up to ~ 
-3(-4) at = 2(3). 

The difference between the (e -|- fe) and fc-corrections 
gives the intrinsic galaxy luminosity evolution, or AMz = 
e-corrcction(2:). This quantity represents the luminosity 
brightening of a galaxy at Ao = \obs x {1 + z)^^ with re- 
spect to its -? = luminosity at the same wavelength. The 
e-corrections at z = 1 and 3, as observed in the B and K 
bands are given in Table 4. The amount of brightening up 
to 2 = 1 is relatively small, as required by the observed 2 
distributions. For 2 ~ 3, AB^ ~ —6 for E/SO galaxies and 
ABz ~ — 1.6 for Sab-Sbc galaxies. In the K band the lu- 
minosity evolution is signiflcantly lower than in the other 
bands (AK.^i ~ -0.6 and AK.^s ~ -2.3 for E/SO galax- 
ies), showing how the concept of mild luminosity evolution 
depends on the wavelength range being considered. 

3.3 Number counts 

Optical counts have reached very deep levels thanks to CCD 
cameras. IR array detectors have made possible to extend 
the faint galaxy work to the K band. Observations in the 
K band are particularly important because light at 2/im is 
relatively insensitive to dust extinction. At high 2 the K 
band samples the well understood rest frame optical range, 
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1 2 3 4 5 6 

z 





K 








Tj-SC 


















// T^Q — SC 






cons — sp 






















(b) 





1 2 3 4 5 6 

Z 



Figure 2. (a) bj band (e + fe)-corrections (solid lines) and 
fc-corrections (dotted lines) for galaxies of different morpholog- 
ical types derived from the BC93 models listed under Q = in 
Table 3. (b) Same as (a) but for the K band. 



Table 5. Slope 7 at faint magnitudes 



Band 


mag. range 




'PLE 


^) 


U 


20-25 


0.49 


0.49 


0.22 


bj 


20-25 


0.45 


0.46 


0.30 


^/ 


20-25 


0.37 


0.36 


0.21 


I 


19-23 


0.34 


0.35 


0.17 


K 


18-23 


0.26-0.30 


0.29 


0.10 



" Prom the n = reference PLE model 

^ Predicted slope at faint limits: 24 < m < 27 



in which spectral evolution is less significant than in the UV 
region sampled at high 2 by the optical bands (Table 4). K 
band counts are thus less sensitive than optical counts to the 
evolution of the stellar population and to the details of the 
SFR and IMF. K counts can therefore be used, at least in 
principle, as a direct test of cosmological models. Because 
of these reasons in this paper we pay special attention to 
fitting at the same time all known properties of faint galaxy 
samples, from the JJ to the K band. 

Figures 3a-e show the differential number counts in the 
{U,bj,rf,I,K) bands. As in most observational papers, 
from which the data have been taken, N{m) is plotted per 



6 Pozzetti, Bruzual and Zamorani 





Figure 3. Differential galaxy number counts per square degree per half magnitude interval as a function of apparent magnitude. The 
sources of the observed data points are indicated in each panel. The lines show the predicted counts for different models: dotted line: 
nE model (Q = 0, = 4.5); solid line: PLE model (Q = 0, Zf = 4.5); dashed line: PLE model (H = 1, zj = 10). (a) U band, (b) bj 
band, (c) rf band, (d) I band, (e) K band. 



half magnitude bin. The sources of the data points are indi- 
cated in the figures. Even though in a few cases the counts 
from different groups show a relatively large scatter, the 
slope 7 of the log N — m relation is reasonably well defined 
in all bands. At bright magnitude 7 ~ 0.6 is close to the 
Euclidean value. Table 5 shows that at fainter magnitude 7 
decreases with increasing filter effective wavelength (Jones 
et al. 1991; Gardner et al. 1993; Djorgovski et al. 1995). 
Despite the large difference in 7 between bj and K, the two 
bands are sampling the same galaxy population. We verified 
that the observed bj counts can be reproduced (in number 
and 7) by just shifting the observed K band counts accord- 
ing to the {B — K)med at a given K given by Gardner et al. 
(1993). The lines in these figures represent the predictions 



of the models that are described in detail below. 
3.3.1 nE model 

The well known excess of faint galaxies above the nE model 
prediction is confirmed (Fig 3). In nE models the SEDs of 

distant galaxies are represented by SEDs that match those 
of nearby galaxies (Table 3), which correspond, on average, 
to old stellar populations. The amount of UV flux produced 
by Post-AGB stars in E/SO galaxies is not as large as that 
attained at early ages. The large fc-correction moves these 
galaxies towards fainter magnitude and, hence, they do not 
contribute to the predicted counts even at the faint limit 
reached by current observations. 
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Figure 4. Differential galaxy number counts as a function of bj 
apparent magnitude. The lines show the predicted counts for two 
different models: solid line: NLE model (Q = I, Zf = 10, r] = 1); 
dashed line: PLE model (0 = 1, 2/ = 10). 

The excess of galaxies with respect to our nE model 
(with f2 = 0) is about a factor of 2 at m 24 in all bands, 
except for the K band in which the nE model is consistent 
with the data. Even at fainter magnitude the observed data 
are never more than a factor of three higher than the nE 
model. The excess reported by other authors (Broadhurst 
et al. 1988; Maddox et al. 1990; GRV90) is larger by a 
factor ^ 2 than our result. This difference is mainly due 
to our choice of normaUzing the model counts to the ob- 
served number in the range 19.0 < bj < 19.5 (§2.3), and to 



our selection of realistic SEDs for early-type galaxies. Our 
SEDs reproduce reasonably well the UV flux of local galax- 
ies, which translates into significant differences between our 
fc-corrections at high z and those used in the quoted papers. 
On the other hand, our results are in reasonable agreement 
with the nE model of KGB93, supporting their conclusion 
that only a moderate amount of spectral evolution is re- 
quired by the data. 

In an n = 1 universe the difference between data 
and model (not shown) at faint magnitude is substantially 
higher, being a factor of ~ 3(4) at bj ~ 24(26). 

3.3.2 PLE models 

From Figures 3a-e we see that the f2 = reference 
PLE model reproduces reasonably well the observed galaxy 
counts in the five bands over a wide magnitude range. The 
most significant discrepancies are with respect to the bright 
bj counts of Maddox ct al. (1990, §2.3), and in K with re- 
spect to the Mobasher et al. (1986) counts and the Gardner 
et al. (1993) HWS counts. At if > 20 the model agrees 
well with the data from Soifer et al. (1994) and McLeod et 
al. (1995), but overestimates the Gardner et al. (1993) HDS 
data and the Djorgovski et al. (1995) counts for K < 23. Ta- 
ble 5 shows the excellent agreement between the observed 
and predicted values of the slope 7 in the five bands. At 
fainter magnitude than the current limits the model predicts 
a significant flattening of N(ni) in all bands (last column of 
Table 5). The faintest r/ and / band counts of Tyson (1988) 
do not show this flattening but show an excess (with large 
error bars) with respect to the model predictions. If faint 
counts keep increasing with a steep slope, then higher values 
oi Zf are likely needed to fit the data. For increasing Zf, the 
fiattening predicted at faint magnitude shifts towards even 
fainter magnitude. Deeper data are necessary to test these 
predictions in detail. 




Figure 5. Contribution of each morphological type to the total differential galaxy number counts as a function of apparent magnitude. 
The lines show the predicted counts for the reference f2 = PLE model. Total counts: solid line; E/SO: dotted line; Sab-Sbc: short-dashed 
line; Scd-Sdm: long-dashed line; vB: dot-dashed line, (a) bj band, (b) K band. 
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Figure 6. bj — rf colour distribution for different bj bins. The solid histograms show the observed distributions from Metcalfe et al. 

(1991, 1995a) and, for the faintest bin, from Koo & Kron (1992; data kindly provided to us by C. Gronwall). When needed, the data 
have been converted to the bj and rf passbands of Couch & Newell (1980). The dotted histograms show the predictions of the f2 = 
reference PLE model. To aid the eye a dashed line is drawn at bj — rf = 0.8. 



In Figures 3are we see that the $7=1 model significantly 
underestimates the faint counts in all bands. Although to 

a lesser extent, this is even true in the K band, where the 
influence of evolution is less important. The deficiency in the 
counts for the n = 1 model is due to the decresising comoving 
volume available for increasing Q. The differences between 
the predictions of this model and the data are minimized 
assuming a high Zf. Even for Zf = 10 the excess in the 
observed counts is ~ 3 at bj = 26. 

3.3.3 Number luminosity evolution (NLE) model 

Models that obey the = 1 constraint of the inflation sce- 
nario require physical effects not considered in simple PLE 
models. Among several possibilities, strong galaxy merger 
events at early ages have been suggested (Rocca-Volmerange 
& Guiderdoni 1990 (hereafter RVG90); Broadhurst et al. 
1992). To test qualitatively this possibility, we have con- 
structed a NLE model in which the LF is allowed to evolve 
with z as proposed by RVG90 

<t>{L, z) = (1 -f z)'>[L(l + z)\ z = 0]. (3) 

This function simulates the merging of faint high z gala^K- 
ies to form bright local galaxies, while conserving the total 



comoving mass density. We have used the same galaxy evo- 
lution models as in the PLE models (Table 3). This simple 
and crude NLE model fits the galaxy counts quite well in all 
bands. A good fit is obtained for 77 = 1 in (3). Fig 4 shows 
the resulting hj counts. More realistic models that simu- 
late the change of the photometric properties of the merg- 
ing galaxies have to be computed (cf. Fritze-v.Alvensleben 
& Gerhard 1994). 

3.3.4 Contribution of each galaxy type to the total counts 

Fig 5 shows the contribution to the total bj and K 
counts of each galaxy type in the f2 = PLE model. In 
6j and U (not shown), early spirals are the dominant popu- 
lation over most of the magnitude range. Evolved E/SOs 
contribute only ~ 20% to the counts at 20 < bj < 22. 
Young E/SO galaxies at their maximum SFR period con- 
tribute ^ 50% and produce the hump in the total counts in 
the range 23.5 < bj < 25.5. This hump has been noticed by 
Metcalfe et al. (1995a) in their counts. E/SO galaxies are 
also dominant for K ^ 22. There is no K hump because 
even a,t z = Zf — 4.5 the K filter samples the blue spectral 
region, not reaching the UV rest frame and missing the sig- 
nature of the high star formation episode. The hump shifts 
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to brighter magnitude for higher Q and to fainter magnitude 
for higher Zf. The onset of the bj hump may seem partic- 
ularly steep in our models because we use a coarse grid of 
galaxy types. With a finer grid of galaxy spectra the bright 
galaxy phases appear more gradually and the onset of the 
hump is loss evident (GK95). The hump in the data (Fig 
3b) suggests some degree of discreteness in the distribution 
of galaxy types, intermediate between ours and GK95's. 

3.4 Colour distributions 

The colour distributions derived from faint galaxy sur- 
veys show a gradual trend towards bluer mean colours at 

fainter magnitude. The result that bj — Vf (photographic) 
becomes significantly bluer beyond bj ~ 22 (Kron 1980), has 
been confirmed by various groups (Shanks ct al. 1984; Tyson 
1984; Infante, Pritchet & Quintana 1986). With CCD detec- 
tors colour distributions have been extended up to bj ~ 27 
(Koo & Kron 1992 and references therein; Metcalfe et al. 
1995a). The median B — K colour of if -selected field galax- 
ies also becomes bluer beyond K ~ 17.5 (Gardner et al. 
1993). To make the comparison of the data and models 
more realistic, we have applied a gaussian error function 
with (J = 0.15 mag to the colour vs. z relations before do- 
riving the colour distributions. This procedure is expected 
to take into account, at least to first order, both the obser- 



vational errors in the colours and the intrinsic dispersion in 
the colours of galaxies of the same morphological type. 

Fig 6 shows the bj —rf colour distribution for different 
bins of bj. At the brightest and faintest bins, panels (a), 
(d)-(f), the agreement between the data and model distri- 
bution is quite good. At intermediate magnitude, panels (b) 
and (c), the model distribution is redder than the data by 
~ A{bj — Vf) = 0.4. Despite this discrepancy, the model 
predicts that galaxies with bj — rf < 0.8 appear at bj ~ 23, 
and are essentially absent at brighter magnitude. 

Fig 7 shows the bj — rf fractional colour distribution 
predicted by the model for the same bj bins of Fig 6 dis- 
criminated by galaxy morphological type. While at bright 
magnitude, panels (a) and (b), the various morphological 
types arc reasonably well separated in colour, at the faintest 
bins, panels (e) and (f), galaxies of all types share the same 
colours. In particular, galaxies bluer than bj — rf = 0.8 are 
not just vB galaxies or late type spirals, but ^ 50% of them 
in panel (c) and about 60% in panel (d), arc high z young 
E/SO galaxies with Zmed ~ 1-5. These are the same galax- 
ies responsible for the hump discussed in §3.3.4. The E/SO 
galaxies are responsible for both the blue (at high z) and 
the red (at low z) tails of the distribution in panels (c) and 
(d). 

Fig 8 compares the Hawaii K-hand survey (Songaila et 
al. 1994) B — K colour distribution with our model. This 
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Figure 8. B — K colour distribution for different K bins. The solid histograms show the observed distributions from the Hawaii K 
band survey (Songaila et al. 1994). The dotted histograms show the predictions of the reference = PLE model. 
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Figure 9. Median colours as a function of apparent magnitude. 
The bottom part of the figure shows (bj — vs. hj. The 

top part shows (B — i<')med '^s. K. The data points are from 
the compilations by Metcalfe et al. (1995a) and Gardner et al. 
(1993). The lines show the predictions for the = reference 
PLE model (solid line), O = 1 reference PLE model (dashed line), 
Q = nE model (dotted line). 

sample has been obtained as a combination of a number 
of different K magnitude limited surveys. The size of the 
surveyed area was diminished as a function of hmiting K 
to provide a roughly constant number of galaxies in each 
magnitude bin from K = Vi Xo K = 20. The model has 
been scaled to the number of galaxies observed in each bin. 
Also in this case the agreement between model and data is 
rather good. 

Pig 9 shows the observed (?>j ?^/)med and (B - K) 
median colours as a function of bj and K. The = PLE 
model reproduces rather well the data. High-a star forming 



galaxies are responsible for blueing the median model colour 
at faint K. This is one possible reason why galaxies with 
colour as red as local E galaxies are not detected at faint 
K (Gardner 1995). The nE mode fails to reproduce the 
data, especially in {B — K)med vs. K, where the difference 
between the PLE and the nE models is large. Even though 
in the = 1 PLE model we have used intrinsically redder 
galaxy SEDs than in the = model (Table 3) , the median 
B — K for the = 1 model is bluer than for the data. The 
main reason for this shift towards bluer colours in the = 1 
model is the younger age of the galaxies in this cosmology, 
despite the higher Zf. 

3.5 Redshift distributions 

Multi-object spectrograph surveys of large samples of faint 
galaxies have reached faint enough magnitude for evolution 
to be detectable. At bj <^ 21 large numbers of z > 1 galax- 
ies are predicted by PLE models in which galaxies undergo 
strong luminosity evolution (Broadhurst et al. 1988, 1992). 
The failure of the z surveys to reveal these galaxies sets 
strong upper limits on the amount of luminosity evolution 
that can be invoked to explain the excess in the number 
counts and the blueing of the median colour discussed in 
the previous sections. Rather than being evidence against 
all PLE models and in favour of the nE model, we argue be- 
low that the observed z distributions of blue-selected faint 
galaxies are consistent with PLE models in which galax- 
ies evolve mildly in luminosity and rule out only the large 
amounts of luminosity evolution assumed in early PLE mod- 
els. 

3.5.1 bj band 

Fig 10 shows the z distributions derived from four surveys. 
Panels (a)-(c) correspond to bj limited samples, while the 
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Figure 10. z distributions for bj selected samples. The source for the observed distributions (solid histogram) is indicated in each 
panel together with Nju = number of galaxies with measured z, N^o-ID = number of galaxies for which no z could be measured, and 

the completeness limit of the survey. The area of the rectangle in each panel = Nno—ID- Dotted line: nE model; solid line: PLE model; 
both for SI = 0, 2/ = 4.5. The model predictions have been scaled to the total number of objects, Af/o + N^^—ID- 



data in panel (d) are a combination of seven fields with dif- 
ferent limiting magnitude, ranging from 23 to 24 in bj, and 
completeness limit > 70% (Glazebrook et al. 1995a). In the 
latter case, our simulation includes the different bj limiting 
magnitude for each field. The data in panel (b) correspond 
to the LDSS survey (Colless et al. 1990), supplemented with 
additional z measurements in a subsample of the original 
LDSS survey areas by Colless et al. (1993). The z distri- 
bution of the f2 = reference PLE model shown in panels 
(a) and (b) is consistent with the data, and does not show 
a significant high z tail. The mean z for the PLE model 
distribution in these two panels are {z) = (0.23, 0.32), in 
excellent agreement with the measured values (0.22, 0.31). 
Colless et al. (1993) reduced the z incompleteness in their 
LDSS subsample to 4.5%. From an analysis of these data 
they conclude that at most ~ 4% of galaxies with bj < 22.5 
have z > 0.7. The PLE model predicts ~ 3% of galaxies 
with 2 > 0.7 in this magnitude range. 

At fainter magnitude, panels (c) and (d), the PLE 
model predicts a high z tail which is not seen in the data. 
Following Glazebrook et al. (1995a), we compare in Table 
6 the values of Zmed (median z) and fo.j (fraction of galax- 
ies with z > 0.7) for their survey and our model. For each 
magnitude range in Table 6, the first line gives Zmed and /0.7 
for the galaxies with measured z, the second line gives up- 
per limits computed by assuming that all the z-unidentified 
galaxies are at ^; > 0.7, and the third line lists the predic- 



Table 6. Glaaebrook et al. (1995) z surveys 



bj range 


source 


^med 


/0.7 


22.5 - 23.5 


data 


0.46 


0.13 




UL« 


< 0.48 


< 0.24 




PLE 


0.50 


0.27 




095* 


0.76-1.31 


0.55-0.82 


22.5 - 24 


data 


0.46 


0.12 




UL 


< 0.56 


< 0.38 




PLE 


0.59 


0.40 




G95 


0.83-1.39 


0.61-0.84 



" Upper limits if 2- unidentified galaxies are at 2 > 0.7 
PLE models by Glazebrook et al. (1995) 

tions of the n = reference PLE model. The model can 
be reconciled with the data only if all, or at least most of 
the galaxies without measured z in the indicated magnitude 
range are at high z. This assumption will be discussed in 
§3.6. Our PLE model is in good agreement with the spectro- 
scopic data of the two _B-sclcctcd SSA13 and SSA22 Hawaii 
Survey fields (Cowie, Hu & Songaila 1995): all 9 galaxies 
with S < 23 have z < 0.7, and for 23 < B < 24, 8 galaxies 
out of 21 (38%) with measured 2 have 2 > 0.7. 

The nE model distribution is in good agreement with 
the data in the four panels of Fig 10. The nE model requires 
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the z distribution of the ^-unidentified galaxies not to be 
very different from that of the galaxies with measured z. 
However, it is important to recall here the failure of the nE 
model to reproduce the counts and the colour distributions 
at faint magnitude (§3.3.1). 

Note that we have not included in our models the 
(1 + dependence of the surface brightness with z, which 
decreases the fraction of expected high z galaxies (Yoshii 
& Peterson 1995). The magnitude of this effect is a func- 
tion not only of the intrinsic parameters of the galaxies, 
but also of the details of the data reduction procedure 
adopted in the construction of the photometric catalogs from 
which the galaxies to be observed spectroscopically are se- 
lected. Other effects, such as dust extinction inside galaxies 
(GK95), which have not been considered in our models can 
decrease the predicted number of high z galaxies. 

To compare different PLE models, we list in Table 6 
the range of Zmed and fo.T computed by Glazebrook et al. 
(1995a) from their PLE models. The Glazebrook et al. val- 
ues are significantly higher than ours. The same is true for 
the GRV90 and Metcalfe et al. (1995a) PLE models. The 
basic reasons for this difference axe the milder galaxy evo- 
lution assumed in our PLE models, and/or the more crude 
treatment of luminosity and spectral evolution of galaxies 
in the quoted papers. For example, in the Glazebrook et 
al. (1995a) model, a single SFR history is adopted for all 
morphological types. On the other hand, the KGB93 and 
GK95 models produce values of Zmed and /0.7 slightly lower 
than ours and, therefore, in good agreement with the data. 
However, the KGB93 model also predicts a significant pop- 
ulation of low luminosity galaxies at 2 < 0.2, which is not 
seen in the Glazebrook et al. data. 

3.5.2 I and K hands 

In Fig 11 we compare the z distribution derived from the 
large and deep z surveys of Crampton et al. (1995, I- 
selected) and Songaila et al. (1994, if -selected) with the 
models. The z distribution in panel (a) is characterized by 
{z) = 0.56, Zmed = 0.57, and fo.r = 0.34. The = ref- 
erence PLE model predicts (z) = 0.58, Zmed = 0.54, and 
fo.7 = 0.30, in excellent agreement with the data. The 
high z E/SO galaxies predicted at bj ^ 23.5 have a mean 
{bj — I) < 1, corresponding to (7) ^ 22.5, and are, there- 
fore, not expected in significant numbers in this survey lim- 
ited at / = 22.0. The nE model z distribution, (z) = 0.48, 
Zmed = 0.47, and /0.7 = 0.18, does not reproduce well the 
tail of high z galaxies observed in this survey. The I band 
LF derived from the Crampton et al. (1995) data set (see 
Fig 6 of Lilly et al. 1995b) , as well as the E galaxy LF con- 
structed from the HST Medium Deep Survey (estimating z 
photorrictrically, Im et al. 1995), show clear evidence of evo- 
lution ioT z < 1. The amount of evolution is consistent with 
the prediction of our PLE model: AI^ ~ —1.2 up to 2 = 1. 

To compare the models with the Songaila et al. (1994) 
z distribution, we have taken into account the decreasing 
area vs. increasing limiting K relationship of the survey as 
discussed in §3.4. The completeness of this spectroscopic 
survey is ~ 100% for K < 18 and ~ 70% for 18 < K < 20. 
At variance with our z distribution for fej -selected samples 
(Table 6), the /0.7 upper limits derived from the AT-selected 
data are significantly lower than the model values. For the 



Crampton et al. (1995) : 




Songaila et al. (1994) : 

K < 20 

N,„=255 




Figure 11. Redshift distributions for / and K selected samples. 
The source for the observed distributions (solid histogram) are 
indicated in each panel. The Crampton et al. (1995) data were 
kindly provided to us by L. Tresse. The two values for complete- 
ness of the Songaila ct al. sample refer to X < 18 and K > 18, 
respectively. Dotted line; nE model; solid line: PLE model; both 
for n = 0, Zf = 4.5. 

observed z distribution in Fig 11, panel (b), /0.7 = 0.08. 
The PLE model predicts /0.7 = 0.23, well above the upper 
limit /0.7 < 0.13 derived assuming that all galaxies with no 
measured z have z > 0.7. This result indicates that there 
may be a problem in one or more of the ingredients or as- 
sumptions that go into our PLE models. The spectroscopic 
incompleteness of the sample is important at the faint level, 
but it cannot be the only cause of the problem since the dis- 
crepancy is already present in the bin 17 < K < 18, where 
the spectroscopic survey is complete. In this bin, the upper 
limit from the observed distribution is /0.7 < 0.16, and the 
model predicts /0.7 = 0.37. We think that the evolutionary 
rate in the K band has to be revised in order to explain 
the observations. We are currently exploring galaxy evo- 
lution models based on different sets of evolutionary tracks 
and stellar spectral libraries (Bruzual 1995) to evaluate their 
effect on PLE models. 

3.6 On the nature of z-unidentified galaxies 

Since most z surveys are based on spectra which cover the 
range from ~3700 to ~7500 A, the [Oil] A3727 fine (most 
prominent feature in low S/N galaxy spectra) can be used 
as a 2 indicator only up to 2 ^ 1. Higher 2 galaxies are 
difficult to identify because of the lack of features in the UV 
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Figure 12. (a) B — K colour as a function of z for the Songaila 

ct al. sample. The solid dots correspond to the ^-unidentified 
galaxies in the survey and are plotted arbitrarily at 2 = 3.5 The 
lines represent the colour predicted for each morphological type 
by the BC93 models listed in Table 3 for the H = reference 
model. E/SO: solid lines; Sab-Sbc: short-dashed line; Scd-Sdm: 
long-dashed line; vB; dot-dashed line, (b) B — I vs. I — K for the 
same sample. The lines have the same meaning as in (a). The 
z = 1 locus is indicated. See text for details. 

rest frame spcetra of typical galaxies. This fact, by itself, 
suggests the possibility that a large fraction of the gala^xies 
with no measured z is at high z. This assumption is also 
supported by an analysis of the colour distributions of these 
galaxies. We discuss here in detail the behavior of the B — I 
and I — K colours in the Songaila et al. (1994) sample. Fig 
12a shows B — K vs. z for all galaxies in the sample, dis- 
criminated by morphological type accordingly to the closest 
proximity of the galaxy colour to the BC93 model colour for 
the class at the same z. Most of the z-unidentified galax;- 
ies (black dots in the plot) are consistent with the model 
colours of E/SOs at z > 0.2. In Fig 12b the same galaxies are 
plotted in the B — I vs. I — K plane. The lines show the ex- 
pected location of galaxies of different morphological type at 
increasing z. Two things arc interesting to note: first, galax- 
ies of all types fall reasonably well in the appropriate region 
in this plane. Second, a large fraction of the z-unidcutified 
galaxies has colours consistent with being carly-typo galax- 
ies at z > 1.0. We have verified a posteriori the reliability of 
the photometric z estimates by examining the measured z of 
galaxies with colours close to the expected colour at « = 1: 



7 out of 8 such galaxies have 0.64 < z < 1.16. Thus, it 
seems reasonable to assume that most of the z-unidentified 
galaxies in this sample are actually at high z. 

This conclusion is consistent with the Keck telescope 
spectroscopic observations of faint galaxies {K < 20, 1 < 
22.5, B < 24.5) by Cowie et al. (1995b). Their data show 
strong evidence of a significant fraction of high-« luminous 
star forming galaxies. Of 333 galaxies that have been ob- 
served, z could be measured for 281 galaxies. 91 galax- 
ies (~ 32%) have z > 0.7 and 40 galaxies (~ 14%) have 
z > 1. Cowie et al. argue that inspection of the B — I vs. 
I — K colour-colour plane suggests that most of the remain- 
ing unidentified objects are luminous high z star forming 
galaxies. 

The same type of analysis cannot be applied to the hj- 
selected survey of Glazebrook et al. (1995a) because only 
B — R\s available for these galaxies. The photometric z esti- 
mates using a single colour are less accurate. The B — R dis- 
tribution suggests a higher percentage of late typo galaxies 
in this sample than in the if -selected sample, in agreement 
with our model for the number counts (§3.3.4). 

4 CONCLUSIONS 

The detailed comparison of a large amount of faint galaxy 
survey data and the predictions of new models has provided 
the following results. 

(i) The n = nE model fits reasonably well the observed 
{U,bj,rf,r) counts up to ~ (21,22,22,23), respectively. At 
fainter magnitude the counts show an excess with respect to 
the nE model in all bands, except in K. The excess in the 



bj counts amounts to Nobs/NnE 



3 at 6j ~ 26, and is lower 



than in previous studies by a factor of ~ 2. This difference 
is due to our choice of normalizing the model counts to the 
observed number in the range 19.0 < bj < 19.5, and to the 
higher UV flux in the BC93 models of local E/SO's than in 
the SEDs used in previous nE models (GRV90, Maddox et 
al. 1990). 

(a) A PLE model, in which galaxies axe characterized 
by the LF and SEDs hsted in Table 1 and Table 3 for the 
(0 = 0, Zf = 4.5, tg = 16 Gyr) Friedmann cosmology, 
provides good fits to most of the existing data from faint 
galaxy surveys. In particular, it reproduces well the faint 
galaxy counts in the U, bj, Vf, I and K bands, as well as 
the colour distributions of blue and Jf-selected samples up 
to the faintest observational limits, and the z distributions 
of bj-, I- and K-selected samples up to bj < 23,/ < 22, 
and K < 17. In the magnitude range 23 < bj < 24, the 
predictions from the model can still be reconciled with the 
data if most of the galaxies without a z determination are 
at relatively high z (z > 0.7). 

(Hi) The steep slope of the APM counts at bright bj re- 
quires strong evolution even at relatively small z (Maddox et 
al. 1990) . This is in conflict with the mild amount of evolu- 
tion required to fit the z and colour distributions at fainter 
magnitude. It has been suggested that the low counts at 
bright magnitude may be produced by either photometric 
errors (Metcalfe et al. 1995b) and/or by the fact that pho- 
tographic surveys may have missed a significant population 
of low surface brightness galaxies because of the relatively 
high isophotal limits in these surveys (McGaugh 1994; Fer- 
guson & McGaugh 1995). 
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(iv) PLE models in a flat Q, = 1 universe cannot reproduce 
several aspects of the data. In particular, the faint counts 
are significantly below the observed ones in all bands, includ- 
ing K. A NLE model in an n = 1 universe, in which galaxy 
density increases as (l + z)'^ due to merger events is found to 
reproduce satisfactorily the counts (see also RVG90). How- 
ever, we consider that a more realistic model, which takes 
into account the change of the galaxy photometric proper- 
ties when the merging occurs, must be explored before this 
result can be taken at face value. 

(v) The Scalo IMF, being relatively poor in massive stars, 
produces models for early-type galaxies which are consid- 
erably less luminous in the UV at early epochs than the 
Salpeter IMF models. This lower flux translates into a 
milder spectral evolution, as required by the observed z dis- 
tributions of faint galaxies. Scalo IMF models are thus pre- 
ferred for early-type galaxies over Salpeter IMF models. 

(vi) The observed bj — Vf and B ~ K colour distributions 
are reproduced more closely by our PLE models if we assume 
that star formation in E/SO galaxies extends over a longer 
period of time (n and T2 models) than in a short box-type 
burst models {B\ model). 

(vii) At faint magnitude {bj <; 23, K ^ 18) our PLE model 
predicts a signiflcant fraction of high-2 galaxies, mostly blue 
young early-type E/SO and Sab-Sbc galaxies. The predic- 
tions of this model are in excellent agreement with the B — K 
colours of the Hawaii Tf-selected survey which failed to re- 
veal galaxies at high z with colours as red as those of local 
E/SO galaxies (Gardner 1995). On the other hand, recent 
results on the LF of red and blue galaxies in the CFRS sam- 
ple suggest evolution of the blue galaxies and not of the red 
galaxies (Lilly et al. 1995b). This appears to be in con- 
tradiction with our model and to favor models in which a 
signiflcant amount of the luminosity evolution needed to fit 
the faint counts is due to spiral rather than early-type galax- 
ies. A discussion of these models can be found in Campos 
& Shanks (1995). 

(via) The most significant remaining discrepancy between 
our PLE model and the data is in the z distribution of the 
Songaila et al. (1994) iC-selected sample &t K > 17. The 
fraction of galaxies at high z predicted by the model is sig- 
nificantly higher, by a factor of ~ 2 — 3, than observed. The 
discrepancy persists even if, as suggested by their colours, 
most of the galaxies which have been observed spectroscop- 
ically and for which no z could be measured, are at high z. 
Thus, on the one hand, the z distribution of these galax- 
ies seems to rule out evolution. On the other hand, the K 
band LF (Glazebrook et al. 1995b) indicates that galaxies 
ai z = 1 are ~ 0.75 magnitude brighter than local ones, 
consistent with PLE models. Only additional spectroscopy 
of faint iiT-selected galaxies can resolve this apparent con- 
tradiction between different data sets. 

In summary, simple PLE models, together with an ap- 
propriate selection of galaxy evolution models (SFR and 
IMF) which provide mild luminosity evolution at least up 
to 2 = 1 in an S7 = universe, can still be considered as 
baseline models. Additional improvements should be con- 
sidered in this framework, for example: (a) Introduction 
of different zj for galaxies of different morphological type 
and/or different luminosity. This could help in explaining 
the colour- luminosity relation, observed in cluster galaxies 



(see Bower, Lucey & Ellis 1992 for ellipticals and Gavazzi 
1993 for both early and late types), (b) Spectrophotomet- 
ric models which include chemical and dynamical evolution 
(Bressan, Chiosi & Fagotto 1994) predict that E/SO galaxies 
approach their local colour at younger ages than in chem- 
ically homogeneous models. This faster rate of evolution 
might bring the predictions of the Q. = 1 PLE model in bet- 
ter agreement with the observations and should be explored 
in detail, (c) Introduction in the models of at least the 
most important observational selection effects, e.g. the sur- 
face brightness selection (McGaugh 1994, Yoshii & Peterson 
1995). (d) A detailed analysis of the effects of extinction by 
dust (GK95, Campos & Shanks 1995), taking into account 
the time evolution of the amount of dust inside galaxies, 
which can be particularly important at high z. 

The interested reader may request the k— and (e -I- 
fc)— corrections required to reproduce these models via e- 
mafl from G.B.A. or L.P. 
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